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Particle number fluctuations are studied in the ideal pion gas approaching Bose-Einstein conden- 
sation. Two different cases are considered: Bose condensation of pions at large charge densities 
Pq and Bose condensation at large total densities of pions p K . Calculations are done in grand 
canonical, canonical and microcanonical ensembles. At high collision energy, in the samples of 
events with a fixed number of all pions, N n , one may observe a prominent signal. When N n in- 
creases the scaled variances for particle number fluctuations of both neutral and charged pions 
increase dramatically in the vicinity of the Bose-Einstein condensation line. As an example, the 
estimates are presented for p + p collisions at the beam energy of 70 GeV. 
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Long time ago in 1924 the Bose statistics was discovered [jl|], and one year later the phe- 
nomenon of Bose-Einstein condensation (BEC) ^ was predicted. Tremendous efforts were re- 
quired however to confirm BEC experimentally. The atomic gases are transformed into a liquid 
or solid before reaching the BEC point. The only way to avoid this is to consider exttemely low 
densities. At these conditions the thermal equilibrium in the atomic gas is reached much faster 
than the chemical equilibrium. The life time of the metastable gas phase is stretched to seconds 
or minutes. This is enough to observe the BEC signatures. Small density leads, however, to small 
temperature of BEC. Only in 1995 two experimental groups succeeded to create the 'genuine' BE 
condensate by using new developments in cooling and trapping techniques [^]. Leaders of these 
two groups, Cornell, Wieman, and Ketterle, won the 2001 Nobel Prize for this achievement. 

Pions are spin-zero mesons. They are the lightest hadrons copiously produced in high energy 
collisions. In the present letter we argue that the pion number fluctuations may give a prominent 
signal of approaching the BEC point. In fact, there is the BEC line in a plane of pion density 
and temperature. The pion system should be in a state of thermal, but not chemical, equilibrium 
to reach the BEC line. This can be achieved by selecting the samples of events with high pion 
multiplicities. Multipion states are formed in high energy nucleus-nucleus collisions, as well as in 
the elementary particle ones. There were several suggestions to search for BEC of 7T-mesons (see, 
e.g., Ref. [Q]). However, complete statistical mechanics calculations of pion number fluctuations 
have never been presented. There is a qualitative difference in properties of the mean multiplicity 
and of the scaled variance of multiplicity fluctuations in different statistical ensembles. The results 
obtained with grand canonical ensemble (GCE), canonical ensemble (CE), and microcanonical 
ensemble (MCE) for the mean multiplicity approach to each other in the large volume limit. This 
reflects the thermodynamic equivalence of the statistical ensembles. Recently it has been found [|5|- 
|9|] that corresponding results for the scaled variance are different in different ensembles, and this 
difference is preserved in the thermodynamic limit. To extract the matter properties from analysis 
of event-by-event fluctuations, one needs to fix the samples of high energy events, and choose the 
corresponding statistical ensemble for their analysis. This is discussed below (see also [^]). 

Let us start with a well known example of non-relativistic ideal Bose gas. The occupation 
numbers, n p , of single quantum states, labelled by 3-momenta p, are equal to n p = 0, 1, .. . In 
the GCE their average values, fluctuations, and correlations are the following [|I0|]: 

(np) = Unl 1 N 7, ((An p ) 2 ) = (n p )(\ + (n p )) = vl, (An p An k ) = v 2 8 pk , (1) 



exp 



1 



where An p = n p — (n p ), m denotes the particle mass, T and fx are the system temperature and 
chemical potential, respectively (throughout the paper we use the units with h = c = k = 1). The 



average number of particles in the GCE reads [ |10| ] : 

V r p 2 dp 



(N) = N{V,T,n) = J> 



p/ " 2tt 2 



exp 
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M )/T 
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where V is the system volume. We consider particles with spin equal to zero, thus the degeneracy 
factor equals 1. In the thermodynamic limit, V — > °°, the sum over momentum states is trans- 
formed into the momentum integral, £ p . . . = (V /2k 2 ) Jq . ..p 2 dp. This substitution, assumed in 
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all formulae below, is valid if the chemical potential in the non-relativistic Bose gas is restricted 
to fx < (or n < m in relativistic formulation). When the temperature T decreases at fixed par- 
ticle number density p = N/V, the chemical potential pt increases and becomes equal to zero at 
T = T c , known as the BEC temperature. At this point from Eq. (||) one finds, N(V,T = T c ,pi = 
0) = l/[mr c /(27T)] 3/2 £(3/2), where £(3/2) ^ 2.612 is the Riemann zeta-function. This gives, 

T c = 2tt[£(3/2)]- 2 / 3 ^— * 3.31 ^— . (3) 
m m 

At /j. = and T < Tc, a macroscopic part, Nc (called the BE condensate), of the total particle 
number occupies the lowest energy level p = 0. At /X = and T <T C the GCE average number of 
particles in the BE condensate is equal to Nc = N[l — (T/Tc) 3 ^ 2 ]. 

Introducing AN = N— (N) one finds the particle number fluctuations in the GCE, 

The limit — fi/T S> 1 gives (n p ) <C 1. This corresponds to the Boltzmann approximation, and then 
from Eqs. (g|) it follows: N(V,T,n) = V exp(n/T)(mT /2k) 3/1 and to = 1. When ju increases 
the scaled variance oo becomes larger, oo > I. This is the well known Bose enhancement effect for 
the particle number fluctuations. From Eq. (Q) at pL — ► one finds oo — > <». Thus, the anomalous 
particle number fluctuations appear in the GCE formulation when the system approaches the BEC 
point. Two comments are appropriate here. First, for finite systems oo remains finite, and oo = °° 
emerges from Eq. (|j) at jx = in the thermodynamic limit V — » °°, when the sums over p are trans- 
formed into the momentum integrals. Second, the anomalous fluctuations of the particle number at 
the BEC point correspond to the GCE description. In the CE and MCE, the number of particles ,/V 
in a non-relativistic system is fixed by definition, thus, £0 c e , = oo mx , e , = 0. 

The average values of the occupation numbers in the relativistic ideal gas of pions equal to: 

M = exp[(7FT^| - Hj)/T] - 1 ' (5) 

where index j enumerates 3 isospin pion states, k + , k~, and 71°, the energy of one-particle states is 
taken as, £ p = (p 2 +m|) l / 2 with m n = 140 MeV being the pion mass (we neglect a small difference 
between the masses of charged and neutral pions). The inequality < m n is a general restriction 
in the relativistic Bose gas, and fij = m n corresponds to the BEC. In Ref. [||] we discussed in 
details the Bose gas with one conserved charge in the CE (V,T,Q = const), i.e. the 7r + 7T _ -gas 
with fixed electric charge. This corresponds to the GCE (V,T,Hq), thus, in Eq. (|J) ju + = \Iq and 
;U_ = — Hq for 7T + and %~ , respectively. Approaching the BEC of n + at \Iq — > m n , one finds the 
relation between Tc and pg = p + — p_ (see Fig. |l|, Left). The picture of BEC of %~ at Q < and 
Hq — > —m n is obtained by a mirror reflection. BEC starts at T = Tc when \Lq = Hq"* = m n . It 
gives: 

Tc tn 2 1 

Pe( T = T cHQ = m x ) = — j^2L- K 2{ nm n/ T c)smii(nm n /Tc), (6) 

71 «=i n 
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Figure 1: Left: The phase diagram of the relativistic ideal Bose gas. The solid line shows Bose condensation 
temperature as a function of the conserved charge density pg given by Eq. (6) at jig = m n . The dashed line 
shows Pq = p+ — p_ = at jlQ = 0. Right: The scaled variances (Q+ e in pion Bose gas, are shown as 
functions of jX* = Hq/T. The solid lines present (0+ e at m* = m n /T = 0.01, 0.1, 0.3, 0.5, 1. The vertical 
dotted lines /x* = m* demonstrate the restriction /I* < m* in the Bose gas. The dashed horizontal line 
presents a value of £(2)/ £(3) ~ 1.368 which is an upper limit for wi~ e reached at fi* = m* — > 0. The 
crosses at jj,* = m* correspond to the points of Bose condensation. The crosses at jj,* — correspond to 
ft)* e (jit* = 0,m* — > 0. The dashed line corresponds to (0^ e in classical (Boltzmann) pion gas (see Ref. [||| 
for details). 



where K2 is the modified Hankel function. At Tc/m n <C 1, Eq. (|h gives, Tq 



In 



K(3/2)p/3 PQ 



2/3 _i 



33\pQ ? 'm n l , which coincides with the non-relativistic formula (|3|). In the ultrarelativistic limit, 
Tc/m K > 1, Eq. (|) gives 



1/2 -1/2 

' 1-1-1 ' 



m 



(J) 



The results presented in Eqs. (^-0) are well known (see, e.g., Ref. [11]). For this system, the 
particle number fluctuations near the BEC line within GCE and CE were recently studied in 
Ref. The scaled variance co + = ((AN + ) 2 ) / (N + ) in the GCE goes to infinity. This is sim- 
ilar to the non-relativistic case. On the other hand, the scaled variance for negative particles, 
ft) = ((AN-) 2 } /(N-), remains finite and even decreases with \Iq. The pion numbers jV + and 
fluctuate in the both GCE and CE. However, the exact conservation imposed in the CE on the 
system charge, Q = N + — AL, suppresses anomalous fluctuations at the BEC point: co^ e (see Fig. [j], 
Right) is finite with the upper limit, £(2)/£(3) ^ 1.368 (see details in Ref. |@]). 



A formation of the pion system with large electric charge density pg in high energy collisions 
does not look realistic. In what follows we discuss a rather different pion system which may be 
created in high multiplicity events [0]. We use the MCE (V,E,Q = 0,N n = const) formulation, 
where the total system energy E, electric charge Q = N + —N- = 0, and total number of pions, 



N X =N +N++ N- 



(8) 
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will be fixed. Such a system can be also described in the GCE (V, T \}Xq=0,}X k ) formulation, with 
H + = \L n + Hq, /i_ = \i n — \Iq, and /lo = P-% in Eq. @. We restrict Hq=0 and consider BEC when 
p. K —nn K . The Hq=0 corresponds to zero electric charge, Q=0 or N+ = 7V_, in the pion system. 

The pion density is equal to p„(T,H n ) =Lpj( n pj) /V ■ The phase diagram of the ideal pion 
gas in p n — T plane is presented in Fig. || (Left). BEC starts at T = Tq when \x n = }i™ ax = «% . It 
gives [0| : 

3Tcm 2 °° 1 

Prc(^ = 7c,jU w = m^) = - 71 V -^2 {nm n /T c ) exp(nm K /T c ). (9) 

2tt 2 

Note an essential difference between Eq. @ and Eq. (Q): a presence of exp(nmx/Tc) in Eq. 
instead of siah(nm % /Tc) in Eq. @. The Eq. @ gives the BEC line shown by the solid line in Fig. || 
(Left). If Tc/m n <. 1, from Eq. (@) one finds, T c = 2^:[3C(3/2)]- 2 / 3 p^ /3 m 3r 1 ^ 1.59p^ /3 /n- *. This 
again corresponds to the non-relativistic limit (||) discussed above, but with a degeneracy factor 
g n = 3. In the ultrarelativistic limit, Tc/m K ^>l, Eq. @ leads to a new relation [0]: 

r c =[7r 2 /3C(3)] i / 3 P y 3 =i.4py 3 , do) 

which differs from Eq. (0) and does not include the dependence on m K . 

Let us consider the region in p n — T plane between the \i n = and \i n = m n lines. The lines of 
fixed energy density, e(r,ju re ) = Lp,y e p ( n p,j)/V, are shown as dotted lines in Fig. || (Left) inside 
this region for three fixed values of e. An increase of p n at constant e leads to the increase of \i % and 
decrease of T. In this letter we discuss how the system approaches the BEC line {pL n = m ]l , T = 7c). 
and do not touch the region (p, n = m K ,T < Tq) below this line where the non-zero BE condensate 
is formed. The GCE (V,T,^Iq,^), MCE (V,E,Q,N K ), and CE (V,T,Q,N K ) are equivalent for 
average quantities, including average particle multiplicities, in the thermodynamic limit. Thus, 
Eq. @ and phase diagram in Fig. 1 remain the same in all statistical ensembles. However, the pion 
number fluctuations are very different in different ensembles. Before starting to calculate the pion 
number fluctuations let us make several comments. 

As an example we consider the high multiplicity events in p + p collisions at IHEP (Protvino) 



accelerator with the beam energy of 70 GeV (see Ref. [12] on the experimental project "Thermal- 
ization", team leader V.A. Nikitin). In the reaction p + p — > p + p + N n with small final proton 
momenta in the c.m.s., the total cm. energy of created pions is E = y/s — 2m p = 9.7 GeV. The trig- 
ger system designed at JINR (Dubna) selects the events with N n > 20 in this reaction. This makes it 
possible to accumulate the samples of events with fixed N n = 30 4- 50 and the full pion identification 
during the next 2 years Jl3|]. Note that for this reaction the kinematic limit is N™ ax = E /m n = 70. 



We stress that the IHEP experiment will measure the both charged and neutral pions (see Ref. []14|]). 
A reliable measurement of A^o number in each event is a crucial point for the identification of the 
BEC line suggested in this letter. The BEC signatures discussed below become useless if the 71° 
number cannot be measured reliably. 

The pion system in the thermal equilibrium is expected to be formed for high multiplicities. 
The volume of the pion gas system is estimated as, V = E / e{T ,p, n ), and the number of pions equals 
to N n = V p n {T ,p, n ). The values of N K at \i n = and \i n = m n for 3 different values of energy 
density e are shown in Fig. @ (Left) for the fixed total pion energy of E = 9.7 GeV. Note that the 
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Figure 2: Left: The phase diagram of the pion gas with /ig = [^||. The dashed line corresponds to 
Pk(T,IJ.k= 0), and the solid line to BEC. The dotted lines show the states with fixed energy densities: 
£ = 6,20,60 MeV/fm 3 . The N n numbers in the fi gure correspond to ji n = and ji n = m n at these energy 
densities for the total pion energy, E = 9.7 GeV. Right: The scaled variance of neutral pions in the MCE is 
presented as the function of the total number of pions [^J. Three solid lines correspond to different energy 
densities: e = 6,20,60 MeV/fm 3 . The total energy of the pion system is assumed to be fixed, E = 9.7 GeV. 
The vertical dotted lines correspond to the points on the BEC line at the specific values of the energy density. 



statistical approach to hadron production in p+p collisions has been used successfully to calculate 



the particle number ratios within the CE [15] and MCE [16]. Such an approach is usually applied to 
a sample of the minimum bias events. Our suggestion has two new points. First, it is a selection of 
the sample of events with high pion multiplicity N n . Most part of the available energy is then spent 
to the pion production. Thus, a strong longitudinal motion seen in the inclusive data is suppressed 
in high multiplicity events because of the energy conservation. The pion system may approach 
the state of global thermal equilibrium with the thermodynamical parameters close to the BEC 
line. To search the BEC effects we suggest to study the specific event-by-event fluctuations of the 
number of pions. This is a second new point of our suggestion. The typical expected temperature 
of the pion gas approaching the BEC line is about of T = 60 — 90 MeV (see Fig. 1). One can then 
calculate the average thermal energy per particle, a rough estimate gives: 3T /2 = 90 — 135 MeV. 
On the other hand, a pion from the p -meson decay has a much larger 'kinetic energy' of about 
245 MeV in the p -meson rest frame. The pion energy becomes even larger due to a presence of 
non-zero rho-meson momenta in the c.m.s. of the p+p collision. Thus, we conclude that for the 
high multiplicity events discussed in this letter a presence of large number of resonances decaying 
into pions is strongly suppressed because of the energy conservation. 

For Q = 0, the average pion multiplicities, (iVo) = {N±) = N K /3, are the same in all statistical 
ensembles for large systems. This thermodynamic equivalence is not, however, valid for the scaled 
variances of pion fluctuations. The system with the fixed electric charge, 2 = 0, the total pion 
number, N K , and total energy of the pion system, E, should be treated in the MCE. The volume 
V is one more (and unknown) MCE parameter. The calculations below are carried out in a large 
volume limit, thus, parameter V does not enter explicitly in the formulae for the scaled variances. 
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The microscopic correlators in the MCE (V,E, Q = 0,Nn) equal to (see also Refs. [0-0]): 



(An J v An{) m . c . e 
where #+=1, q-= — l, qo=0, 



\A\ 



qjqiMgg+Mmc + s pi s k jM £e - (e pi + e kj )M ne 



pj 



<«p> (1 + <«p» 



(n p ) = {exp[(Jp> + ml-n s )/T] - 1}" 



(11) 



(12) 



|A| is the determinant and My are the minors, 

M qci = A(K 2 )A(e 2 )-(A(7ie)) 2 , M m = A(q 2 ) A(e 2 ), M ee = A(q 2 ) A{k 2 ), M ne = A{q 2 ) A(jre), 

(13) 

of the correlation matrix A, 

fA(q 2 ) 
A = A(7T 2 ) A(ne) 
^ A(ks) A(£ 2 ) 

The matrix A ( |I4| ) has the following elements: 

T 2\ 



(14) 



A(^ 2 )=I^ p 2 J = 2£ Up 2 , 



A(£ 2 ) = ^g 2 ^. = 3£^Uj 



2„,2 
P 



P,J 



A(^)=IUp 2 y - = 3£ Up 2 , 
pj p 

pj p 



(15) 



Note that the first term in the r.h.s. of Eq. (|llj) corresponds to the GCE. Correlations between 
differently charged pions, j^i, and between different single modes, p/k, are absent in the GCE: 



CO 



CO 



CO 



CO 



1 + 



Ip(" P ) 2 

E p ("p) 



(16) 



similar to non-relativistic result but with (n p ) given by the relativistic relation (fl2|). In the 
GCE the numbers A^ + , Af_, and A^o fluctuate independently of each other. The Bose effects in 
the pion system are small if \i n = 0. For \i % = 0, one finds CO = 1.0 1 1. 12 in the temperature 
interval T = 40— 160 MeV (note that CO = 1 in the Boltzmann approximation). The Bose effects 
increase with \i n , and CO — > °° at \i n — > m n , i.e. approaching the BEC line the GCE calculations 
give anomalous fluctuations for A^ + , A^_, and A^o- 

The MCE (V,E, Q = 0,N K = const) formulation means the restrictions of the exactly fixed total 
system energy E, electric charge Q = N + —N- = 0, and total number of pions N n ^ for the each 



microscopic state of the system. This changes the pion number fluctuations. From Eq. ( |1 1| ) one 
notices that the MCE fluctuations of each mode p are reduced, and the (anti)correlations between 
different modes p ^ k and between different charge states appear. This results in a suppression of 
all scaled variances COm. c . e . in comparison with the corresponding ones CO in the GCE. A nice feature 



of the MCE microscopic correlators (1 1 ) is that although being different from that in the GCE, they 
are expressed with the quantities calculated in the GCE. The MCE scaled variances depend on two 
GCE parameters: T and \i K . 
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The substitution of Eqs. (12-15) in Eq. ( JTTj ) and straightforward calculations lead to the fol- 
lowing MCE scaled variance for neutral pions [0] : 

llp,k( An p An k) m.c.e. 2 

<,,. = = 3 « - ( 17 ) 



where ft) is given by Eq. (16) and corresponds to pion fluctuations in the GCE. Due to the condi- 
tions, N + = N- and N+ +N- +N = N n , and equal average multiplicities, (No) = (N+) = (N-) = 
Nx/3, it follows [0]: 

®m.c.e. = 4 <c* = g W > and C = J ^-e- = 3 W ' (18) 

where A^/, = N + +N-. The behavior of (O mce (17) is shown in Fig. || (Right). To make a corre- 
spondence with Njt values, we consider again the p + p — > p + p+N n collisions at the beam energy 
of 70 GeV and take the pion system energy to be equal to E = 9.1 GeV. Despite of the MCE sup- 
pression the scaled variances for the number fluctuations of 71° and jt^ increase dramatically and 
abruptly when the system approaches the BEC line. 

The following inequalities are always hold for particle number fluctuations in different en- 
sembles: 0nh.ce. < G>c.e, < (Og.c.e.- Therefore, if the anomalous BEC fluctuations are present in the 
MCE, they are also exist (and even larger) in the CE and GCE. The reverse statement is not true. 
The anomalous BEC fluctuations of the GCE may disappear in the CE or MCE. We found that 
for the system with N n = const and Q = the anomalous BEC fluctuations are not washed out by 
exact conservation laws of the CE and MCE. This is an advantage of the system with N n = const 
and 2 = 0. Let us repeat again (see discussion just after Eq. (0) and Ref. [§]) that the anoma- 
lous BEC fluctuations at high charge density pg disappear in the CE and/or MCE. As another 
instructive example let us consider the MCE (V,E ,Q=0,N c h = const), i.e. fixed N c h = N + +N-, 
instead of fixed N n ([|). The corresponding GCE formulation gives the following pion chemical 
potential: /i + = \i n , Hq = 0, = \i n in Eq. (||) (/ig = 0, as before, because of Q = condition). 
When /ijr — * m % the system approaches the BEC line for 7i + and n~ . The thermodynamic behavior 
and position of this BEC line can be easily found. Approaching the BEC line one can also find 
ft)* — > 00 in the GCE. The pion number fluctuations are, however, very different in the both CE 
(V, T,Q = 0,N c h) and MCE (V,E, Q = 0,N ch ). In the statistical ensembles with fixed N ch and Q no 
anomalous BEC fluctuations are possible. The numbers of and jV_ are completely fixed by the 
conditions Q = N + - = and N c h = N + + AL = const, thus, (S)f e = C0^ c e = 0. The number 
No fluctuates, but jj, = 0, thus, neutral pions are far away from the BEC line and their fluctuations 
are small, ft) 1, in all statistical ensemble formulations. 

The broad distributions over No and N c h close to the BEC line also implies large fluctuations 



of the / = No /N c h ratio. These large fluctuations were suggested (see, e.g., Ref. Qi7|]) as a possible 
signal for the disoriented chiral condensate (DCC). The DCC leads to the distribution of / in the 
form, dW(f)/df = l/(2y/J). The thermal Bose gas corresponds to the /-distribution centered 
at / = 1/2. Therefore, /-distributions from BEC and DCC are very different, and this gives a 
possibility to distinguish between these two phenomena. 

The calculations presented in this letter should be improved by taking into account the finite 
size effects, pion-pion interactions, and some other effects. Examples discussed in Refs. [Q, ||] 
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demonstrate that the thermodynamical limit for the average multiplicities and scaled variances 
is reached rather quickly. Normally, if the average pion multiplicity is about of N n = 10, the 
deviations from thermodynamic limit in the ideal pion gas are only a few percents. The Bose- 
Einstein condensation is a phase transition phenomenon. Thus, the infinite volume limit is of a 
principal importance. A strict mathematical meaning of the phase transition (and its order) has only 
sense in the infinite volume limit. Of course, the real systems are finite. For our applications this 
means that the scaled variance for neutral pions shown in Fig. 2 does not increase up to 'infinity', 
but it is restricted from above in the finite system. A detailed study of the finite size effects for 
the BEC is now under investigation and will be published elsewhere. The BEC temperature of 
about T = 60 — 90 MeV corresponds to the pion number density of p n = 0.1 — 0.15 fm" 3 (see 
Fig. 1). This particle density is not too large (it is smaller than the normal nuclear density). This 
probably may justify the ideal pion gas approximation considered in the present letter as a first 
step in the modelling of multi-pion states. The effects of pion interactions will be discussed in 
the future studies. Preliminary estimates suggest that the BEC signatures suggested in the present 
letter may survive the complications. A crucial point is the analysis of the samples of high N K 
events. The required N K values for the BEC are much larger than the average pion multiplicity per 
collision, thus, these high N n events are rather rare and give negligible contributions to inclusive 
observables in high energy collisions. With increasing of N % in the sample with fixed total energy, 
the temperature of the pion system has to decrease and it approaches the BEC line. This can happen 
in different ways: at constant energy density e, at constant pion density p n , or with decreasing of 
both £ and p n . The pion system should move to the BEC line one way or another. In the vicinity 
of the BEC line (no BE condensate is yet formed) one observes an abrupt and anomalous increase 
of the scaled variances of neutral and charged pion number fluctuations. This could (may be even 
should) be checked experimentally. 
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